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ABSTRACT 


To facilitate the study of engagement scenarios with high energy lasers, the 
Directed Energy Physics Group at the Naval Postgraduate School developed a laser 
performance code called Atmospheric NPS Code for HEL Optical pRopagation 
(ANCHOR). This code uses well-known analytical scaling laws and a scriptable user 
interface to allow the quick exploration of multi-dimensional parameter studies. Recently, 
a new capability was added that incorporates relative platform / target motion. This study 
demonstrates this new capability, and compares ANCHOR results with those obtained 
with the full diffraction code built into WaveTrain. 
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I. 


INTRODUCTION 


High Energy Laser (HEL) technology has evolved to the point that laser weapons 
now have the capacity to be used in naval platforms. Their advantages over conventional 
weapons in military applications include unlimited magazines (as long as they are fed 
with power), low cost per shot, and very precise delivery of the payload at the speed of 
light. 

However, unlike conventional weapons, laser weapon performance is greatly 
affected by weather conditions and atmospheric properties. Thus, the atmospheric 
propagation of HELs needs to be studied in order to characterize their effectiveness and 
performance as laser weapons. 

In this thesis, the effects of relative target and platform motion on the atmospheric 
propagation of HELs is studied, in particular how this relative motion affects an 
atmospheric phenomenon called thermal blooming. This study is conducted by modifying 
a scaling code called ANCHOR that was developed by the Directed Energy Physics 
Group at the Naval Postgraduate School to include motions of the platform and the target. 
We begin this study with an overview of directed energy weapons in Chapter II. Next, we 
discuss the atmospheric propagation of high energy lasers in Chapter III. In Chapter IV, 
we discuss how the performance of the laser weapons is characterized and how the 
relative motion affects the atmospheric propagation. Then, we study different cases in 
order to investigate how the propagation is affected by the relative motion in Chapter V. 
Einally, we conclude our study Chapter VI. 
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II. OVERVIEW OE DIRECTED ENERGY WEAPONS 


A. HISTORY 

Directed energy (DE) weapons are used in order to damage a target by depositing 
energy onto it for a specific dwell time. The most studied and promising version of 
directed energy weapons is lasers (light amplification by stimulated emission of 
radiation). 

The idea of a laser started with Albert Einstein in 1916 as he predicted how 
stimulated emission at particular wavelengths can be achieved . The first laser, the ruby 
laser, was demonstrated by Ted H. Maiman in May 1960, which had around 1 watt of 
output power [1]. After this accomplishment, the importance of this new technology was 
realized and further research led to the first gas laser (He-Ne laser) demonstration in 
December 1960, the demonstration of the first solid-state laser (Nd:YAG) in 1964, the 
demonstration of the first carbon dioxide laser in 1964, and the demonstration of the first 
free electron laser in 1977 [1]. 

The technology in this area evolved to the point that laser weapons now have the 
capacity to be used in naval platforms. A solid state laser (SSL) weapon was tested in 
April 2011 on the USS Paul Foster (DD 964) in the Pacific Ocean, tracking and setting 
fire to multiple, small unmanned boat targets. The laser system was used at high power 
mode more than 35 times and withstood actual maritime conditions including eight-foot 
waves, winds of 25 knots, and fog [2]. 

B. ADVANTAGES 

In order to understand the advantages of directed energy weapons, we need to 
look for the reason why they are wanted to be used in military applications and what are 
their capabilities that conventional weapons do not have. 

The DE weapons apply damage at light speed with very good precision. The 
speed can overcome the target’s ability to maneuver kinetic weapons, and the precision 
enables damage to a specific part of the target that can enhance the lethality while also 
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reducing the chance of collateral damage. The precision is at the level that the main beam 
of the Boeing Airborne Laser, a megawatt-class chemical laser, was able to hit targets 
that were 500 kilometers away with pinpoint accuracy [3]. 

Another important advantage of electric DE weapons is an essentially unlimited 
magazine if the platform can support the weapon with power, unlike a conventional 
weapon that has a specific amount of munition. This is especially important in swarm 
attacks where conventional weapons may be exhausted before all the targets are 
eliminated. Furthermore, the cost per shot for a DE weapon is cheaper than a 
conventional weapon. “The shipboard fuel needed to generate the electricity for firing an 
electrically powered laser would cost about a dollar per shot” [4]. Finally, the effects of 
the engagement on a target can be varied by changing the power of the DE weapon. As a 
result of this, nonlethal and lethal applications may be supported by a single weapon [1]. 

C. TECHNOLOGIES 

1. Solid State Lasers 

An SSL uses a solid substrate that contains dopants as a gain medium. These 
dopants are excited by a pumping source in order to achieve a lasing transition. A cavity 
that consists of two mirrors, one reflective and one partially transmittive, is used to 
increase the intensity of the light produced by the lasing transition via the process of 
stimulated emission. 

Two kinds of SSLs that will potentially be used in military applications in the 
near term are fiber lasers and slab lasers. Slab lasers have a gain medium formed as a 
rectangular slab, while a fiber laser uses a long cylindrical medium consisting of an inner 
and an outer cladding. Fiber lasers have the advantage of easier cooling than slab lasers 
due to their geometry, while more light coupling can be achieved with slab lasers due to 
their large aspect ratio. 

Since they are a more mature technology, and relatively compact compared with 

Free Electron Eases (FEEs), SSEs are planned to be used in DE weapon applications in 

the near term. However, unlike FELs, they are restricted to operation at a specific 

wavelength, which changes with the material used as the gain medium. Also, due to the 
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heat produced inside the solid gain medium operation, the possibility of damaging the 
medium exists and limits the output power levels of SSLs. FELs do not have this 
problem. Also, in order to reach the desired DE application levels for SSEs, multiple 
lasers need to be combined, which results in a decreased beam quality and/or increased 
system complexity. 

2. Free Electron Lasers 

Eree Electron Easers (EEEs) are used to harvest the energy of unbounded 
relativistic electrons as coherent light. There are two main types of EEEs used: oscillator 
EEEs and amplifier EEEs. Either type has several core components. An injector and 
accelerator are the first parts, where the electron beam is created and accelerated to 
relativistic speed. Then these electrons pass through the undulator, where a periodically 
changing magnetic field forces the electrons to wiggle and emit light. In an oscillator 
EEE, the light is stored in a resonator and builds up over many passes. In an amplifier 
EEE, the light from a seed laser source is amplified over a single pass. While EEEs are 
inefficient at extracting energy from the relativistic electrons on a single pass, the leftover 
energy retained in the electron beam can be very efficiently recovered to enhance the 
overall wall plug efficiency [1]. 

Since an EEE does not have a solid or gaseous gain medium, thermal effects are 
not as limiting a factor like other types of lasers for both beam quality and output power 
levels. The “basic architecture of an EEE offers a clear potential for scaling up to power 
levels of one or more megawatts” [4]. Oscillator EEEs are not restricted to be operated at 
certain wavelengths, so the wavelength can be tuned to match different atmospheric 
transmission sweet spots. The wavelength of an amplifier EEE is fixed by the seed laser 
used. 

EEEs are typically large, expensive, and the technology used is not as mature as 
SSEs [1]. Another drawback of EEEs is that they create harmful radiation during 
operation that needs to be shielded for both personnel and material safety. This may limit 
the use of EEEs in military applications in near term, but their capability of reaching MW 
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power levels makes them a viable choice for the applications requiring these power 
levels. 
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III. ATMOSPHERIC PROPAGATION OF 
HIGH ENERGY LASERS 


A. ATMOSPHERIC EXTINCTION 

Light loses some of its energy along the path it transverses due to atmospheric 
extinction, which includes scattering and absorption. 

Scattering occurs when photons interact elastically with particles in the 
atmosphere that results in deflection of some of the photons in arbitrary directions. “The 
particles responsible for scattering cover the sizes from gas molecules (~10^*cm) to large 
raindrops and hail particles (~1 cm)” [5]. Absorption is an inelastic process in which 
some of the photons’ energy is transferred to particles, and ultimately results in heating 
the atmosphere. This heating, in turn, can distort the laser beam due to thermal blooming. 
Thermal blooming will be discussed further in Chapter III, Section C. 

The effects of scattering and absorption are quantified by Beer’s Law, 

P{z) = P,e^\ (1) 

In Equation (I), P{z) is the laser power after traveling a distance z , P^ represents the 

initial power at z = 0 (i.e., living the beam director) and s represents the total extinction 
coefficient (in units of inverse length), which consists of absorption and scattering 
coefficients due to molecular and aerosol effects, 

8 = aa+am+ Pm- ( 2 ) 

In Equation (2), a and /? represent absorption and scattering coefficients, while the a 
and m subscripts stand for aerosol and molecular contributions, respectively. Each 
coefficient in Equation (2) “depends on the wavelength of the laser radiation” [6]. In 
order to understand absorption and scattering thoroughly, aerosol and molecular effects 
need to be discussed separately. 

I. Molecular Effects 

Since the wavelengths used in HEEs are in the infrared region, so all of the effects 
need to be studied in this region specifically. 
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The atmosphere is composed of a group of gases with nearly fixed concentration 
and a group of gases with variable concentration [5]. Nitrogen (N 2 ) and oxygen (O 2 ) 

compose ~ 99 % of the atmosphere. However, they are not responsible for the absorption 
in near-infrared region where SSLs operate, since they are homonuclear molecules and 
lack dipole moments. As seen in Figure 1, water vapor (H 2 O) and carbon dioxide (CO 2 ) 

are predominantly responsible for the atmospheric absorption in the infrared region. This 
dominance is the result of having many vibrational modes. Molecular absorption is 
therefore a strong function of the wavelength. 


0, 


H,0 COt H|0 

M|0 CO| 0| 

,111 11 11 1 


CO, 0, H,0 CO, CO, 

I Ml 



Figure 1. Atmospheric transmittance measured over 1820-m horizontal path at 

sea level. Source: [6]. 

When the wavelength of the laser is considerably larger than the particle’s sizes, 
the scattering process is described by Rayleigh scattering. The Rayleigh scattering cross 
section, which is proportional to the likelihood of a scattering event occurring, can be 
approximated as [5] 
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2 


( 3 ) 



where a is the diameter of the molecules and assumed to be much smaller than the 
wavelength of the laser beam (a <sc 2), and n is the refractive index of the molecules. As 
wavelength increases, the effect of Rayleigh scattering on atmospheric propagation 
decreases rapidly (cr^ oc . The sky’s blue color during daytime can be explained by 

Rayleigh scattering; since blue light’s wavelength is shorter than red light’s, blue light is 
scattered more than red light. 


Using the Rayleigh scattering cross section, the molecular scattering coefficient 
can be calculated as [5] 

/^JZ) = aM)]N(r)dr, (4) 

where A(r) is the number density of molecules with radius r. 


2. Aerosol Effects 

Aerosol extinction coefficients are less dependent on the light wavelength than 
molecular extinction coefficients. These coefficients are determined by the combined 
effects of the aerosol’s index of refraction and the size distribution of aerosol particles, 
according to the theory of Mie scattering [7]. In general, the aerosol effects are more 
prominent at shorter wavelengths, at least in the IR spectral region [7]. 

The effect of aerosol scattering on atmospheric propagation is generally far 
greater than molecular (Rayleigh) scattering except at very short wavelength values due 
to molecular scattering’s dependence on wavelength as 2Figure 2 shows that for 
wavelengths up to ~ 0.3 pm, molecular (Rayleigh) scattering dominates the total 
scattering extinction, but aerosol scattering starts to become dominant at near-IR and 
longer wavelengths. 
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Figure 2. Calculated aerosol absorption and extinction coefficients for 23-km 
visibility conditions and a continental aerosol model. Source: [8]. 


B. ATMOSPHERIC TURBULENCE 

Atmospheric turbulence is the result of temperature differences, convection, wind 
shear, and inertial cascades of larger scale turbulences to smaller scales [1]. 
Inhomogeneities in the atmosphere (eddies) result in density variations and, in turn, index 
of refraction differences along the laser beam path. “The shimmering of distant objects on 
a hot day or the twinkling of the stars at night” are examples caused by index of 
refraction variances due to atmospheric turbulence [9]. These index differences also 
distort the wave fronts of an HEL, potentially causing the laser beam to break up into 
beamlets that travel in slightly different directions. 

An important parameter that is used to characterize the severity of atmospheric 
turbulence is the refractive index structure coefficient with units of m^^^^. This 
parameter increases as the temperature differences of the eddies increases, so an increase 
in cl means an increase in atmospheric turbulence. Cl tends to decrease with altitude. 

Typical Cl values near the surface vary from about 10for relatively weak 
turbulence to 10 m for relatively strong turbulence. Furthermore, turbulence varies 
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throughout the day and tends to be smallest “about one or two hours before sunrise and 
after sunset,” while the peak values tend to be measured around noon on sunny days [9]. 
The irradiance patterns for a laser propagating for 5 km through a region of constant 
are shown in Figures 3 through 5. In Figure 3, the turbulence is extremely weak at 
= 10, which has almost no effect on the irradiance. 



X dimension [m] 

2 _ 10 _ 2/3 

Figure 3. Irradiance with weak turbulence (C„ =10 m ' ). Source: [10]. 

The irradiance pattern with stronger turbulence (C^ =10is shown in 

Figure 4. The figure shows that the laser beam gets distorted as the atmospheric 
turbulence increases. 
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Figure 4. Irradiance with moderate turbulenee (= 10). 

Source: [10]. 


The irradiance pattern with even stronger turbulence (C^ =5x10 m is 
shown in Figure 5. The figure shows that the distortion level reaches to the point where 
the peak irradiance is down by over an order of magnitude relative to Figure 3. 
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Figure 5. Fradiance with strong turbulence (= 5 x 10m^^^^). 

Source: [10]. 
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The Fried parameter is the diameter over which transverse beam coherence is 
maintained. It incorporates information about along the path, the wavelength of the 
laser and the propagation distance. Larger values of correspond to weaker turbulence 
and more complete constructive interference at the target. For a focused beam 
propagating along a path of constant , the Fried parameter is given by [1] 


Tq = 0.33 


/I'’/' 




3/5 


( 5 ) 


where A is the wavelength of the laser light and R is the range of the target. Larger 
values and longer ranges result in smaller revalues, while longer wavelengths result in 
longer revalues. Generally, turbulence becomes significant when is smaller than the 
beam director size. As an example, the Fried parameter is «12.5 cm for A = l pm, 
cl = 10and 7? = 5 km. 


C. THERMAL BLOOMING 

Absorption of laser light by the molecules and aerosols in the atmosphere results 
in heating the air, which then changes the density of air along the laser path. These 
density gradients create refractive index changes in this path. These alterations in the 
optical properties of the path cause the air to act like a concave lens, which in turn results 
in defocusing and spreading of the laser beam. This nonlinear atmospheric phenomena is 
called thermal blooming. The effect of thermal blooming is especially crucial in low 
altitude engagements (due to higher water vapor and aerosol content) such as ground/sea 
surface-based scenarios, as well as slow moving targets in light winds due to reduced 
beam slewing. Wind across the beam path tends to reduce the severity of thermal 
blooming since it replaces the air warmed by the laser with cooler air. 

There are two limiting cases of thermal blooming: conduction dominated thermal 
blooming and convection dominated thermal blooming [11]. Conduction-dominated 
thermal blooming exists for the cases where there is basically no wind (v=0). This case 
is rarely met for the HEL applications along the entire beam path, though there can be 
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places along the path that are stagnant. Convection-dominated thermal blooming exists 
when there is the wind or beam motion. This case is met for most HEL applications and 
has been studied most extensively [11]. 

When wind/beam motion exists, cooler air enters from one side of the beam and 
warms as it transverses the beam. This causes the air to heat up in the downwind 
direction, and thus the refractive index drops in the downwind direction. Since light 
bends toward higher refractive indices, the irradiance bends into the wind, often forming 
a distinct crescent shape as illustrated in Figure 6. 




Figure 6. Typical thermal blooming effected irradiance pattern with wind 
coming from left side {P = 100 kW , R = 5 km, = 10 m/s) . 

Source: [10]. 
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IV. METHODOLOGY 


ANCHOR (Atmospheric NPS Code for HEL Optical pRopagation), a scaling 
code developed by the NPS DE group, is used to characterize laser performance. The 
code incorporates parameters like turbulence, wind, jitter, and extinction. Eor this thesis, 
the effect of relative target and platform motion on atmospheric propagation was also 
added to the code. Eor different laser power and laser wavelength inputs, time averaged 
irradiance and power in the bucket plots can be acquired using this code. The code takes 
atmospheric parameters from other models like EEEDR, NAVSEAM, MODTRAN, and 
CO AMPS. 

In order to estimate the irradiance at the target, a couple of parameters need to be 
taken into account. Atmospheric extinction, beam quality, platform jitter, atmospheric 
turbulence, and thermal blooming are all effects that cause the irradiance to decrease 
along the propagation path. Assuming a Gaussian beam at the source, the time-averaged 
irradiance profile, with platform jitter and turbulence-induced wander and spreading 
effects included, is estimated as [7] 

( 6 ) 

where is the total power leaves the beam director, s is the total atmospheric 
extinction coefficient, I is the range to the target, is the time-averaged 1/e radius of 
the mean irradiance profile in the absence of thermal blooming, Sj^ represents the Strehl 
ratio that accounts for the effect of the thermal blooming on the time averaged irradiance. 
The time-averaged beam radius on the target is given by 



In Equation (7), the j , d and t subscripts represent contributions from beam jitter, 
diffraction and atmospheric turbulence, respectively. 
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A. LINEAR EFFECTS ON SPOT SIZE 
I. Beam Jitter 

Beam jitter effect is caused by vibrations and/or tracking errors on the beam 
director platform. The time averaged mean square radial displacement of the focal spot 
due to beam jitter in Equation (7) is given approximately as [7] 

~ ( 8 ) 
where 6^^ represents the angular variance due to jitter and i is the range to the target. 
Typical values of angular jitter are on the order of a few jurad . 


2 . 


Diffraction 


The effect of diffraction on the laser beam spot size for an infinite Gaussian 
source can be estimated as [7] 




k Wn 


+ w„ 


A 


1 - — 

V Fj 


(9) 


where i represents the range to the target, k = Itt/X is the wave number, is the He 

beam radius at the beam director, F is the focal range of the laser beam, and 
represents the beam quality factor. The value is the measurement of how many times 
an ideal beam can fit across the real beam in each transverse direction at the focus. For an 
ideal Gaussian beam, M^=l whereas >1 for non-ideal laser beams [12]. For this 
thesis, the focal range F is equal to the target range £. For this condition, the second 
term in the square root drops out and Equation (9) simplifies to 

M£ 


= 


kw„ 


( 10 ) 


3. Turbulence 

The effect of the atmospheric turbulence (see Chapter III, Section B) on the laser 
beam spot size can be estimated as [13] 


11 

w, «—,for r^<D, 
kr^ 


( 11 ) 
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where is the Fried parameter (see Equation (5)) and D represents the diameter of the 
beam director. The effect of atmospheric turbulence increases with the range of the target 
and decreases with the Fried parameter since it represents the diameter that the 
transverse coherence of the laser beam is maintained. 


B. THE EFFECT OF THERMAL BLOOMING ON SPOT SIZE 


The effect of the thermal blooming on the laser beam propagation through the 
atmosphere is discussed in Chapter III, Section 3. The effect on the laser beam spot size 
can be estimated mathematically using the thermal blooming distortion number, which is 
calculated as [1] 

fMilIhKhffc, (12) 

p path ^wind (z)D(z) ’ 

where z is the distance along the beam path, P is the initial power of the laser, is the 


ambient air density, is the specific heat at constant pressure, a is the absorption 
coefficient, T is the transmission of the laser energy, n^. =dnldT =-(n()-l)/r is the 
change in the refractive index of air with respect to temperature where italic T represents 
the temperature, represents the effective wind speed perpendicular to the laser beam 


and D is the beam diameter. Absorption near the beam director causes more significant 
thermal blooming distortion than the absorption near the target [1]. This effect is called 
lever arm effect. In order to compensate this effect, the following weighting function can 
be included in the integral in Equation (12) [1], 


F{z) = - 


1 + 


T25^ 


V 


J 


z/i 

1-zA 


where D represents the primary aperture diameter in this case. 


(13) 


If there is no change in the atmospheric conditions along the path. Equation (12) 
can be simplified as [1] 

-AyflkPriT-aT . 


No=- 


PoCpDV„,^ 
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(14) 



Thermal blooming effects become significant when >25 [1]; the laser beam gets 
more distorted as increases. In Equation (14), it is seen that the transverse effective 
wind decreases the distortion level due to thermal blooming, and in turn can increase the 
irradiance at the target. The irradiance patterns for a laser (P = 100 kW ) propagated 5 
km through a region with different effective wind speeds are shown in Figures 7 and 8. In 
Figure 7, the effect of the thermal blooming on irradiance is shown for =3 m / s . 

Irradiance xio* 

12 

10 

8 

CM 

E 

6 ^ 

4 

2 

0 

- 0.4 - 0.2 0 0.2 0.4 

X dimension [m] 

Figure 7. Thermal blooming effected irradiance patterns with = 3 m / s 

=465,P = 100 kW, P = 5 km). Source: [10]. 

In Figure 8, the effect of the thermal blooming on irradiance is shown for 
'^wind =5 m/s . It is seen that as the transverse wind speed increases, the irradiance also 

increases (note the different irradiance scales on the graphs), although the crescent¬ 
shaped distortion of the beam is enhanced. 
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Figure 8. Thermal blooming effected irradiance patterns with = 5 m / s 
{N^=219 ,P = \00 kW, i? = 5 km). Source: [10]. 



The Strehl ratio Sj.^, the ratio of the peak irradiance of the thermal blooming 
distorted beam to the peak irradiance of a reference beam with no thermal blooming can 
be calculated as [1] 


5 ^ 

™ h l + 

where a = 0.0625 and m = 2 for a Gaussian distributed beam, whereas a 
m = 1.2 for a uniform distributed laser beam. 


(15) 

0.01 and 


Thermal blooming also affects HEL laser systems such that after a certain laser 
output power, , the irradiance on the target starts to decrease. This critical power 

level can vary with the laser wavelength and the atmospheric and weather conditions, and 
can be estimated by [1] 


P = 

crit 


f p ^ 

r -1 / 

UJ 



1/m 


(16) 
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In ANCHOR, the distortion number from Equation (14) is calculated as well as a 
weighted distortion number which uses the same equation but includes a weighting 
function granting more effect of thermal blooming to the midpoint along the beam path. 
Using these two numbers, an analytic expression for the Strehl ratio is used that is 
calibrated against full diffraction codes. 

1. Effective Transverse Wind Speed 

In order to calculate the effective transverse wind speed across the beam, both 
target and platform motions and the natural wind need to be taken into consideration. As 
a result of the relative target and platform motion, the beam director must slew to keep 
the laser on the target. This slew motion also introduces an effective wind across the 
beam. 

In Figure 9, the parameters needed to calculate the transverse wind speed is given. 
The calculations are made relative to the ship using spherical coordinates. 



The target velocity relative to the ship is calculated as 


V =V -V, 

target target ship ’ 


(17) 


where represents the target velocity and is the velocity of the ship. Likewise, 


the wind velocity relative to the ship is calculated as 
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( 18 ) 


^wind ^wind ^ship ’ 


where represents the natural wind velocity. Angular velocity of the target about the 
ship is calculated using the relative target velocity as 


CO = 


r xV 

target target 


r 

target 


(19) 


where represents the displacement vector of the target relative to the ship. Slew 
velocity along path r is calculated using the angular velocity of the target as 


Vsie.=^xr. 


( 20 ) 


The effective wind velocity is calculated using the slew velocity and relative wind 
velocity as 

yeJf=Kind-ysle.- ( 21 ) 

Finally, the effective transverse wind (i.e., the cross wind), is calculated using the 
effective velocity as 




=(%4)h{%-e)0. 


where (j) and 6 are the spherical unit vectors, so 

= - sin (l)x + co^(l)% 
and 


( 22 ) 


(23) 


0 = cos 0 cos x+ cos 0 sin (jdy - sin Oz. (24) 

We created a MATLAB script using these equations that plots the effective 
transverse wind speed for a constant target range and varying orientation about a circle 
centered on the ship with the same and parameters. In Figure 10, the 

effective transverse wind speed is plotted for =10(m/s)y and 


ywind =-10(m/s)y . In these figures, arrow length represents the magnitude of the cross 

wind along the beam. Also, the color of the arrows evolves from blue to yellow to red as 
the magnitude increases. In this example, there is no cross wind along the vertical axis, 
because in that case the wind is in the same direction as the ship and target motion, so a 
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stagnation zone exists in this direction. This region is where thermal blooming will be the 
most problematic. Along the horizontal axis, the effect of the thermal blooming on 
propagation is the least because the strongest cross wind exists along this path. 



Figure 10. Effective transverse wind speed for = 10(m / s) y and 

^v.,w=-10(m/s)^ 


In Figure 11, the effective transverse wind speed is plotted for = 10(m / s) y, 

^target = 10(tTi / s) X uud = 10(m / s)x + 10(m / s)y . Here, we also see that a stagnation 

zone exists for certain target/ship orientations. Furthermore, the position along the path 
where this stagnation zone occurs changes with the orientation. 
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Figure 11. Effective transverse wind speed for = 10(m / s) y, 
V,arge, = / s) X and = 10(m / s)x + 10(m / s)y. 


C. OTHER PARAMETERS NEEDED TO CHARACTERIZE LASER 

PERFORMANCE 

Other than the irradiance, laser performance can be characterized by power in the 
bucket and dwell time. 


1. Power-in-the-Bucket 


The power delivered to a circle, or “bucket,” of specified radius on the target is 
called power-in-the-bucket. The calculation of the power in the bucket changes with the 
chosen bucket radius r^in comparison with the thermal blooming affected beam radius. 

The thermal blooming affected beam radius can be calculated using the Strehl ratio as 
[14] 




w, 


fot.eff 


'TB 


(25) 


Since ^ accounts for the thermal blooming effect, the irradiance on the target can be 


rewritten as [14] 
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(26) 


(I)= - 

If , power in the bucket can be calculated approximately as [14] 

( 27 ) 

as all of the delivered power falls within the defined bucket. If , then the power 

in the bucket can be calculated approximately as [14] 

PB^{I)[7rr^). (28) 

In this case, we are not interested in all the power delivered to the target. The total 
power delivered to the target is more than the power in the bucket. 

2. Dwell Time 

A laser weapon is used to transfer energy to the target, and in turn damage it by 
melting the material used on the target. The time needed for melting through the target is 
called dwell time. 


a. Power Loss Mechanisms 


In order to calculate the dwell time, power loss on the target needs to be known. 
There are two power loss mechanisms: conductive and radiative losses. These losses 
decrease the power level transferred to the target. 


Conduction loss exists as a result of the heat increase at a specific area on the 
target and the conduction of the heat through the other parts of the target. Conduction loss 
can be quantified using Fourier’s heat conduction law as 


P = 

cond 


I^Kond -Tq) 


Ax 


(29) 


where k is the thermal conductivity of the material, is the melting temperature of the 


material, is the ambient temperature of the material, is the cross-sectional 

surface area of the surrounding material and Ax is the distance of the temperature 
gradient in material. This distance depends on the material properties and the dwell time; 
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for typical metals and dwell times of r^„^;/~10s, A;c~2cm [15]. Typically, the 
conductive power loss is much greater than the radiative power loss. 

Radiative power loss is due to blackbody radiation as the material heats up and 
glows. Using Stefan-Boltzmann law, the power radiated away is proportional to the 
fourth power of the blackbody’s temperature azT^) , the radiated power due to 
temperature change after laser dwell can be calculated as 

= Aea{T:-T:), (30) 

where A is the area over which the radiation occurs, ^is the material emissivity and cr 
is the Stefan-Boltzmann constant (~ 5.670x 10 * W m'^ K'"*). 


Total power loss due to conduction and radiation is calculated as 


p = P + P 

loss cond rad. 


(31) 


The power available to melt the target can be quantified using the power loss as 

p — p X f — P 1321 

melt trans J target ^ loss, ^ * 


where represents the power transferred to the target and is the fractional 

target absorption. For aluminum, is ~ 0.2. In order to melt the target, the power 

transferred to the target must exceed the power loss on the target, otherwise the target can 
not be melted. 


b. Energy Needed to Melt the Target 

The energy needed to melt a target varies with the melting temperature of the 
target material, heat of fusion of the target material, and the mass of the material. The 
energy needed to reach the melting point of the material is calculated as 

= CpinAT, (33) 

where Cp represents the specific heat capacity of the material, m is the mass of the 
material needed to melt and AT =T^-T^^ is the change in the temperature. The energy 
needed to melt the material at the melting point is calculated as 

Q^=mNH, (34) 
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where A// represents the heat of fusion of the target material. The total energy to melt 
the target material is calculated as 

Qn,eU=Q\+Q2- ( 35 ) 


c. Dwell Time Calculation 

Dwell time can be calculated using these calculations as 

_ _ Qmeh 

dwell p 

melt 


(36) 


Dwell time can vary with the power level transferred to the target, the power loss on the 
target, material of the target, and the volume needed to melt. 


d. Laser Performance Calculation Example 

As an example, required on-target irradiances for melting a 100 cm^ circular area 
on an aluminum target are calculated for dwell times of 3 and 6 seconds and material 
thicknesses of 1 and 3 mm. Table 1 shows the physical properties of aluminum used in 
these calculations. 


Table 1. Physical Properties of Aluminum 


Parameter 

Value 

Density 

2.70 

cm 

Thermal conductivity, k 

W 

237 

mK 

Distance of the temperature gradient, Ax 

~ 2 cm 

Melting temperature, 

933.47 K 

Material emissivity, s 

0.05 

Specific heat capacity of the material, Cp 

897 ^ 

kgK 

Heat of fusion of the material. Mi 

400 H 

kg 
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Calculated total energies to melt the target and total power losses for 1 mm and 3 
mm thicknesses are shown in Table 2. 


Table 2. Total Power Loss and Total Energy to Melt Calculation Results 


Parameter 

For 1 mm 
thickness 

For 3 mm 
thickness 

Total energy to melt the target, 

26.14kJ 

78.42 kJ 

Total power loss, 

2.68 kW 

8.04 kJ 


Using the parameters in Table 2, the required on-target irradiances are shown in 
Table 3. 


Table 3. Required On-Target Irradiances to Melt the 
Target Calculation Results 


Parameter 

For 1 mm 
thickness 

For 3 mm 
thickness 

Required on-target irradiance to melt 
target for = 3 seconds 

MW 

1.13 , 

m 

3,4 T 

m 

Required on-target irradiance to melt 
target for = 6 seconds 

0.7 “'f 
m 

. , MW 

2.1 — 
m 
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V. RESULTS 


A. NOTIONAL LASER/ATMOSPHERIC PARAMETERS 

In order to evaluate the effects of relative target and platform motion on the 
atmospheric propagation of HELs, the parameters used as inputs to ANCHOR are given 
in Table 4. 


Table 4. The Input Parameters Used in ANCHOR. 


Laser Output Power P 

200 - 1000 kW 

Laser Wavelength X 

1.064 pm 

Beam Quality 

3 

Platform Height Hp 

10 m Above Sea Level 

Target Height Hp 

0.4 - 1300 m Above Sea Level 

Diameter of Beam at Source D 

0.3 m 

Beam Shape at Source 

Uniform 

Beam Jitter 6^^ 

5x10rad 

Speed of Target 

0 ^ 198 m/s 

Speed of Platform 1/ 

0^10 m/s 

Speed of Wind 

0^3 m/s 

Bucket Radius 

0.05 m 

Fractional Target Absorption 

0.15 
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In order to evaluate the extinction coefficients for X = 1.064 |jm along the laser 
beam path, MODTRAN is used. MODTRAN is a commonly used radiative transfer 
model developed by Spectral Sciences Inc. and the U.S. Air Force Research Laboratory 
[16]. The parameters used as input for MODTRAN are given in Table 5. 


Table 5. Tropical Maritime Input Parameters Used for MODTRAN. 


Parameter 

MODTRAN Model 

Geographical-seasonal 
Model Atmosphere 

Tropical atmosphere (15 ° 
Latitude) 

Aerosol Model 

Maritime extinction 

Cloud and Rain Model 

No rain 

Meteorological Range 

23 km 

Temperature, 

300 K at sea level (varies with 
altitude) 

Absolute Humidity 

19.7 g/m^ (80% relative humidity) 
at sea level (varies with sea level) 


The refractive index structure coefficient along the laser beam path is 
calculated using the Hufnagel-Valley model. This model characterizes the variations of 
cl with the altitude for a given surface value of C,J and is specified by the following: [1] 


C'(/z) = 5.94x10 


-53 


IT 

V27y 


10 -A/IOOO 


1 lu 

n e 


-16 -^1500 


+ Ae 


-/i/lOO 


(37) 


where h is the height in meters, W is the high altitude wind speed in m/s and A is the 
cl surface value. In this thesis, A is set to the moderate value of 10m^^^^. 
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B. COMPARISON WITH WAVETRAIN 


1. Wavetrain 

In order to validate the ANCHOR code results, Strehl ratio outputs are compared 
with the output of WaveTrain for the same conditions. WaveTrain is a software tool 
developed by MZA Associates Corporation “for high fidelity modeling of advanced 
optical systems such as laser weapons systems” [17]. It is a full diffraction code, which 
potentially makes it more accurate than ANCHOR, but ANCHOR is much faster than 
WaveTrain. All of the WaveTrain outputs for this thesis were provided by Dr. Conor 
Pogue. 

In WaveTrain, the Strehl ratios are calculated using the same absorption 
coefficients, target, platform and wind velocities, target and platform heights, and target 
distances with ANCHOR. Strehl ratio values are the division of the peak irradiance with 
thermal blooming and the peak irradiance without thermal blooming for WaveTrain. 

2. Results and Discussion 

In Figures 12 and 13, the comparison of Strehl ratios of ANCHOR and Wavetrain 
is plotted in polar coordinates for a constant target range and varying orientation with the 

same and parameters. Atmospheric turbulence effects are ignored for 

both codes. In Figure 12, the parameters used are = 10(m/ s) y, = 10(m/ s) y, 

=-10(m/s)y , /’ = 100kW,/l = 1.064 pm (typical wavelength used in SSLs), 
R = 5 km and height of the target hj. = 10 m. The effective transverse wind speed plot for 

the same platform, target and wind velocities was shown in Figure 10. It is seen that, 
ANCHOR tends to underestimate the Strehl ratio more as the cross wind along the laser 
beam increases. 
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Strehl Ratio Comparison 
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Figure 12. Strehl ratio output eomparison of Wavetrain and ANCHOR for 
A, = 10(m / s) y , V',,,,, = 10(m / s) y , V.,„ = -10(m / s)y , /> = 100 
kW, X = 1.064 pm , R = 5 km and hj = 10 m. Source: [10]. 

In Figure 13, the parameters used are =10(m/s)y, =10(m/s)x, 

=10(m/s).r + 10(m/s)5', P = 100kW, A =1.064 pm, R = 5km and height of the 
target hj, = 10 m. The effective transverse wind speed plot for the same platform, target 

and wind velocities was shown in Figure 11 . It is seen that the Strehl ratios of the two 
methods agree generally except the regions where ANCHOR slightly underestimates the 
Strehl ratio in the case where the cross wind speed is larger at the zone closer to the 
platform and slightly overestimates the Strehl ratio where the cross wind is larger at the 
zone closer to the target. 
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Figure 13. Strehl ratio output comparison of Wavetrain and ANCHOR for 
y,up = 10(m/ s) y , = 10(m/s)x , = 10(m/ s)x + 10(m/ s)y , 

P = 100 kW, Z = 1.064 pm , R = 5 km and = 10 m. Source: [10]. 


C. ANCHOR EXAMPLES 

The irradiance of the laser will be plotted for a constant target altitude, hj , for 
varying target ranges from the platform at the origin. Every point except the origin 
corresponds to a target with the same velocity. The velocity of the wind also 
does not change for any target orientation. In Figure 14, these target orientations 

are displayed for several points on the polar plot with the velocities of =10(m/s)5', 
y,ar,e, = 10(m/ s)X and = 10(m /s)x + 10(m/s) y. 
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Figure 14. Visualization of target orientation for = 10(m / s) y, 
V.arser = 10 ( iTi / s ) X and =10(m/s)x + 10(m/s)3). 


1. Simple Cases 

In all of the irradiance plots, latitudinal range is along the y-axis, whereas 
longitudinal range is along the x-axis. The contour lines on the plot connect points of 
constant irradiance (outer contour line is for 1 MW/m^ , middle contour line is for 
5 MW/m^ and inner contour line is for 10 MW/m^ ) . The magnitude of the on-target 
irradiance is defined with the color bar scale. In Figure 15, the on-target irradiance is 
plotted for = -3(m /s)y, hj =10m and P = 200kW. There is no relative platform 

and target motion for this plot. Since we have a stagnation zone on the vertical axis, and 
the cross wind gets stronger along the horizontal axis, we have an on-target irradiance 
pattern that is compressed along the vertical axis and elongated along the horizontal axis. 
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-5 0 5 

Longitudinal Range (km) 

Figure 15. Simple case scenario for = -3(m / s) hj =10 m, P = 200 kW 

and /I = 1.064 |jm. 

In Figure 16, the on-target irradiance is plotted for = 10(m/s) 
Kin, =-3(m / s)y, hj = 10 m and P = 200 kW. With the introduction of platform motion, 

the magnitude of the cross wind increases at every point except the stagnation zone for 
this case. The on-target irradiance also increases at every point except the stagnation zone 
(along the vertical axis in this case). 


Figure 16. 



Longitudinal Range (km) 

Simple case scenario for = 10(m/ s) y, = -3(m/ s) y, 
hj = \0m, P = 200kWand/l = 1.064 pm. 
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In Figure 17, the on-target irradianee is plotted for = 10(m/s) 

K..«.r=10(m/s)i), y,,„, =-3(m/s)3), /v=10m and P = 200kW. With the 

introduction of target motion in Figure 17, the cross wind gets stronger at every point 
except the stagnation zone, and in turn the irradianee increases everywhere except along 
this zone. 



-5 0 5 

Longitudinal Range (km) 

Figure 17. Simple case scenario for = 10(m / s) y, = 10(m/ s) y, 

Kind = -3(m/ s) y, /v = 10 m, F = 200kW and A = 1.064 pm. 

In Figure 18, the on-target irradianee is plotted for =10(m/s)y, 

Krget = 55 (m / s) y, Km = -3 (m / s) y, K =l0m and P = 200kW. Since we have more 

target speed in this case, we also have a stronger cross wind than other simple case 
scenarios. As a result of this, we have the largest irradianee at every point except the 
stagnation zone. Since there is no change in cross wind along the stagnation zone, the 
irradianee at this zone stays the same for every simple case scenario. 
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Figure 18. Simple case scenario for = 10(m / s) % Vtarget = 55(m / s) y. 

Kind = -3(m/ s) }), /v = 10 m, F = 200kW and A = 1.064 pm. 

2. Complex Cases 

The on-target irradiances are plotted for more complex cases in this section. In 
Figure 19, the on-target irradiance is plotted for = 2(m/s).r-2(m/s) y, hj=l0m 

and P = 200 kW. There is no relative platform and target motion for this plot. Since the 
wind blows from 315° azimuth, we have a stagnation zone on the line that connects 
315° and 135° azimuth, whereas the cross wind gets stronger through the line that 
connects 45° and 225°. As a result of these zones, we have a on-target irradiance pattern 
that is compressed along the 135°-315° line and elongated along the 45°-225° line. 
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Figure 19. Complex case scenario for = 2(m / s) x- 2(m/ s) /v = 10 m, 

P = 200 kW and A = 1.064 . 

In Figure 20, the on-target irradiance is plotted for =10(m/s)5', 

=2(m/s)x-2(m/s)5', hj=l0m and F = 200kW. With the introduction of 

relative platform and target motion, the magnitude of the effective cross wind changes at 
every point for this case. This results in a new orientation of the on-target irradiance 
pattern. 
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Figure 20. Complex case scenario for = 10(m / s) y, 

=2(m/s)x-2(m/s)y, = 10 m and P = 200kW. 
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In Figure 21, the on-target irradianee is plotted for = 10(m/s) 

=10(m/s)x, =2(m/s)x-2(m/s)j), =10m and P = 200kW. With the 

introduction of target motion in Figure 21, the effective cross wind changes at every point 
again, and in turn this results in a new on-target irradianee orientation. 
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Figure 21. Complex case scenario for =10(m/s)5', y,arget =10(m/s)x, 

K.W =2(m/s)x-2(m/s)y, /v =10m, P = 200kWand 2 = 1.064 

pm. 

In Figure 22, the on-target irradianee is plotted for 1?^;^ =10(m/s)5', 
=55(m/s)x, =2(m/s)Jc-2(m/s)y,/v =10m and P = 200kW. With the 

increase in target speed, both on-target irradianee magnitude and pattern orientation 
changes for every target position on the plot. 
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Figure 22. Complex case scenario for =10(m/s)5', y,arget = 55(m/s)x, 
=2(m/s)x-2(m/s)5), /v=10m, P = 200kW and A = 1.064 

pm. 

D. PRACTICAL EXAMPLES 

In order to characterize the effect of relative target and platform motion on the 
atmospheric propagation of HELs, practical examples will be evaluated via the changes 
in on-target irradiance for different cases in this section. 

1. UAV Target 

In this sub-section, the on-target irradiances will be evaluated with varying UAV 
velocities, varying platform (ship) velocities, varying wind velocities, different UAV 
altitudes and initial laser powers. 

a. The Effect of Target Altitude and Target Speed 

In Figures 23 through 25, the on-target irradiances are plotted for constant laser 

output power, wind velocity and platform velocity which are P = 200 kW, 

^wind = -3(m / s)x and = 2(m ! s)x + 2(m / S)y, respectively. Target velocities 

change from 20(m/s) y to 198(m/s) y . In Figure 23, the on-target irradiances are plotted 

for target height Hj = 300 m. The amount that the atmospheric propagation is effected 

by the relative target and platform motion decreases with increasing speed of the target. 
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As the target speed increases, the effective cross wind starts to be dominated by the target 
speed, and cooling of the laser beam path leads to much less thermal blooming. 
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Figure 23. On-target irradiances for varying target velocities. 

Kind = -3(m / s)V, = 2(m / s)V + 2(m / s) y, P = 200 kW and 

=300m. 


In Figure 24, the on-target irradiances are plotted for Hj^ = 800 m. The same 

trends are seen as in Figure 23. Also, as the target altitude increases, the slew rate of the 
beam director following the target decreases and in turn the effect of relative target 
motion decreases. 
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Figure 24. On- target irradiances for varying target velocities, 

y«in 4 = -3(m / s)x, = 2(m / s)x + 2(m / s)y, F = 200 kW and 

=800m. 


In Figure 25, the on-target irradiances are plotted for =1300 m. Again, the 
same trends are observed as in the previous figures. 
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Figure 25. On-target irradiances for varying target velocities, 

y«in 4 = -3(m / s)x, = 2(m / s)x + 2(m / s)y, F = 200 kW and 

= 1300 m. 


b. The Effect of Laser Output Power 

In Figure 24, the on-target irradiances are plotted for constant target heights, wind 
velocities, and platform velocities, which are Hj =%Q0m, =-3(m/s)v and 

Kft/p = 3(m/s)v-i-2(m/s)y , respectively. Target velocities and laser output powers 

change from 20(m/s).r to 60(m/s)x and 200 kW to 1 MW, respectively. It is seen that 
the effect of the relative target and platform motion on the atmospheric propagation 
depends on the output laser power, since target/platform motion mainly affects thermal 
blooming which is more prevalent at higher laser powers. 
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Figure 26. On-target irradiances for varying target velocities and laser 

output powers, = -3(m / s)x , = 2(m / s)x + 2(m/ s)}' and 

=800m. 


c. The Effect of Target Direction 

In Figure 27, the on-target irradiances are plotted for constant laser output power, 
wind velocity, platform velocity, and target speed, which are P = 200kW, 

Kind =-3(m/s)v,y^;„^ =2(m/s)v-L2(m/s)y and =122 m/s, respectively. Target 

direction and target altitudes change from 0° to 90° azimuth and 300 meters to 
1300 meters, respectively. Figure 27 shows that the direction of the target affects the 
atmospheric propagation such that the orientation of the on-target irradiance contours 
changes in a similar manner with target heading regardless of the target height. 
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Figure 27. On-target irradiances for varying target directions and heights, 
K.W = -3(m / s)v , y, = 2(m / s)x + 2(m / s) y, P = 200 kW and 


^ target = 122 m/S. 


d. The Effect of Platform Velocity 

The effect of the platform velocity on the effect of relative target and platform 
motion on the atmospheric propagation will be studied with different platform speeds and 
platform headings. 
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(1) The Effect of Platform Direction 

In Figure 28, the on-target irradiances are plotted for constant laser output 
power, wind velocity and target height, which are P = 200 kW, = -3(m ! s)x , and 
Hj =500 m, respectively. Target speed changes from 20 m/s to 180 m/s in the y- 

direction and platform direction changes but with a constant speed of 9 m/s. As the target 
speed increases, it starts to dominate the irradiance pattern and the platform direction has 
less of an impact. 



-5 0 5 

Longitudinal Range (km) 



-5 0 5 

Longitudinal Range (km) 



-5 0 5 

Longitudinal Range (km) 





=20(m/s)y, V .. =9(m/s)x 


ro 

^ 0 
(0 
c 

^ -5 

3 ^ 

ro 

-> -10 


Longitudinal Range (km) 

_ =100(m/s)y, V^|^.p=9(m/s)y 


-5 0 5 

Longitudinal Range (km) 

V.ar^et=''80(m/s)9^^,Ws)9 


- I u 


-5 0 5 

Longitudinal Range (km) 



-5 0 5 

Longitudinal Range (km) 



-5 0 5 

Longitudinal Range (km) 



-5 0 5 

Longitudinal Range (km) 


Figure 28. On-target irradiances for varying target speeds and platform 
headings, = -3(m / s)v , P = 200 kW and Hj = 500 m. 
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(2) The Effect of Platform Speed 

In Figure 29, the on-target irradiances are plotted for constant laser output power, 
wind velocity and target height, which are P = 200kW, =-3(rr\/s)x, and 

Hj =500 m, respectively. Target speed changes from 20 m/s to 180 m/s in the y- 

direction and platform speed changes from 5 m/s to 15 m/s in the x-direction. As the 
target speed increases, it starts to dominate the irradiance pattern and the platform speed 
has less of an impact. 
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Figure 29. On-target irradiances for varying target and platform speeds, 
'^wind = “3(m ! s)x, P = 200 kW and Hj. = 500 m. 
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e. The Effect of the Wind Velocity 

The effect of the wind velocity on the effect of relative target and platform motion 
on the atmospheric propagation will be studied with different wind speeds and wind 
headings. 

(1) The Effect of Wind Direction 

In Figure 30, the on-target irradiances are plotted for constant laser output power, 
platform velocity and target height, which are P = 200kW, =-9(m/s)v, and 

Hj =800 m, respectively. Target speed changes from 20 m/s to 180 m/s in the y- 
direction and wind direction changes but with a constant speed of 4 m/s. As the target 
speed increases, it starts to dominate the irradiance pattern and the wind direction has less 
of an impact on the irradiance like the effect of platform direction. 
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Figure 30. On-target irradiances for varying wind directions and target speeds, 

= -9(m /s)x, P = 200kW and = 800 m. 


(2) The Effect of Wind Speed 

In Figure 31, the on-target irradiances are plotted for constant laser output power, 
platform velocity and target height, which are P = 200kW, =-9(m/s)v, and 

=800 m, respectively. Target speed changes from 20 m/s to 180 m/s in the y- 

direction and wind direction changes from 2 m/s to 8 m/s. As the target speed increases, it 
starts to dominate the irradiance pattern and the wind speed has less of an impact on the 
irradiance like the effect of platform speed. 
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Figure 31. On-target irradiances for varying wind and target speeds, 
= -9(m /s)x, P = 200kW and = 800 m. 


2. Small Boat Target 

In this sub-section, the on-target irradiances will be evaluated with varying target 
(small boat) velocities, varying platform (ship) velocities, varying wind velocities, 
varying laser output powers and constant boat height {Hj = 0.4m). 

a. The Effect of Target Velocity 

The effect of relative target and platform motion on the atmospheric propagation 
will be studied with different target speeds and target headings. 
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(1) The Effect of Target Direction 

In Figure 32, the on-target irradiances are plotted for constant wind velocity, 
platform velocity and target height, laser output power, and target speed, which are 

^w,w=-3(m/s)v, y,,,,=5(m/s)y, H,=QAm, P = 200 kW, and m/s, 

respectively. The direction of the target changes for each plot. It is seen that the 
atmospheric propagation is affected in each plot with the change in target direction. 
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Figure 32. On-target irradiances for varying target directions, 

k,w = -3(nn / s)i, V,„„ = 5(m / s)^, ^> = 200 kW, = 9 m/s and 

Hj = 0.4 m. 


(2) The Effect of Target Speed 

In Figure 33, the on-target irradiances are plotted for constant laser output power, 
wind velocity, platform velocity and target height, which are P = 200 kW, 

^wind =-3(m/s)v, =5(m/s)5' and Hj =0.4m, respectively. Target speed changes 

from 0 m/s to 18 m/s in the x-direction. The on-target irradiance plot changes noticeably 
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until the target speed dominates the relative motion. After this speed, the ehange is very 
small. 
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Figure 33. On-target irradiances for varying target velocities, 

= -3(m / S)i , = 5(m / s)y, P = 200kW and H, = 0.4 m. 


b. The Effect of Laser Output Power 

In Figure 34, the on-target irradiances are plotted for constant wind velocity, platform 
velocity and target height, which are =-3(m/s)v, =l{ml s)x + l{xr\ls)y and 
Hj =0.4m, respectively. Target speed changes from 0 m/s to 18 m/s in the x-direction 

and laser output power changes from 200 kW to 1 MW. It is seen that the effect of the 
relative target and platform motion on the atmospheric propagation increases with the 
increase in the output laser power. 
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Figure 34. On-target irradiances for varying target velocities and laser output 
powers, = -3(m / s)x , = 2(m / s)x + (m/s)y and Hj = 0.4 m. 


c. The Effect of Platform Velocity 

The effect of the platform velocity on the effect of relative target and platform 
motion on the atmospheric propagation will be studied with different platform speeds and 
platform headings. 

(1) The Effect of Platform Direction 

In Figure 35, the on-target irradiances are plotted for constant laser output power, 
wind velocity and target height, which are P = 200 kW, = -3(m / s)x , and Hj = 0.4 
m, respectively. Target speed changes from 7 m/s to 18 m/s in the y-direction and 
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platform direction changes but with a constant speed of 9 m/s. Unlike the UAV case, the 
target speed does not dominate the irradiance pattern since the difference between target 
and platform speed is not so large. As a result of this, the on-target irradiance pattern gets 
effected with any change in the platform heading, although this effect decreases as the 
target speed increases. 
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Figure 35. On-target irradiances for varying target speeds and platform 
headings, =-3(m/s)i, P = 200kW and =0.4 m. 


(2) The Effect of Platform Speed 

In Figure 36, the on-target irradiances are plotted for constant laser output power, 
wind velocity and target height, which are P = 200 kW, = -3(m / s)x , and Hj = 0.4 


54 









































































m, respectively. Target speed changes from 7 m/s to 18 m/s in the y-direction and 
platform speed changes from 5 m/s to 15 m/s. Unlike the UAV case, the target speed does 
not dominate the irradiance pattern since the difference between target and platform 
speed is not so large. As a result of this, the on-target irradiance pattern gets affected with 
any change in the platform speed, although this effect decreases as the target speed 
increases. 
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Figure 36. On-target irradiances for varying target and platform speeds, 
V^ind = -3(m / s)x , P = 200 kW and = 0.4 m. 
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d. The Effect of the Wind Velocity 

The effect of the wind velocity on the atmospheric propagation will be studied 
with different wind speeds and headings. 

(1) The Effect of Wind Direction 

In Figure 37, the on-target irradiances are plotted for constant laser output power, 
platform velocity and target height, which are P = 200 kW, = -3(m / s)x , and 
Hj =0.4 m, respectively. Target speed changes from 7 m/s to 18 m/s in the y-direction 

and wind direction changes but with a constant speed of 7 m/s. Unlike the UAV case, the 
target speed does not dominate the irradiance pattern since the difference between target 
and wind speed is not so large. As a result of this, the on-target irradiance pattern gets 
affected with any change in the wind heading, although this effect decreases as the target 
speed increases. 
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Figure 37. On-target irradianees for varying target speeds and wind headings, 

= -3(m Is)x, P = 200kW and = 0.4 m. 


(2) The Effeet of Wind Speed 

In Figure 38, the on-target irradianees are plotted for constant laser output power, 
platform velocity and target height, which are P = 200kW, =-3(m/s).v, and 

=0.4 m, respectively. Target speed changes from 7 m/s to 18 m/s in the y-direction 

and wind speed changes from 0 m/s to 10 m/s. Unlike the UAV case, the target speed 
does not dominate the irradiance pattern since the difference between target and wind 
speed is not so large. As a result of this, the on-target irradiance pattern gets affected with 
any change in the wind speed, although this effect decreases as the target speed increases. 
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Figure 38. On-target irradiances for varying target and wind speeds, 
= -3(m /s)x, P = 200kW and = 0.4 m. 
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VI. CONCLUSIONS 


According to the results, the effects of relative target and platform motion on the 
atmospheric propagation of HELs depend on parameters like target altitude, target 
velocity, platform velocity, laser output power, and wind velocity. As the target altitude 
increases, the effect of the relative motion on the propagation decreases, because the slew 
rate of the beam director decreases with the increase in target altitude. The wind velocity 
can either decrease or increase the relative motion effect on propagation for a specific 
target acquisition, depending on the specific velocities of the platform and target as well 
as their relative positions. As an example, as the target speed becomes faster than the 
platform speed, the relative motion becomes dominated by the target speed. Changes in 
both target and platform direction affect the propagation of the laser regardless of the 
target altitude, target and platform speeds, and the wind velocity, because the effective 
cross wind changes with the change in directions for each case. The effect of the relative 
target and platform motion on the atmospheric propagation also depends on the output 
laser power, since target/platform motion mainly affects thermal blooming, which is 
more prevalent at higher output laser powers. 

As a result, the relative target and platform motion significantly affects the 
atmospheric propagation of HELs, and the level that the propagation is affected depends 
on the parameters discussed above. HELs can be used more effectively by the users 
considering these effects. 
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